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• Abstract 

Ph, We calculate the transverse momentum distribution for the production of 

D . massive lepton-pairs in longitudinally polarized proton-proton reactions at 

collider energies within the context of perturbative quantum chromo dynamics. 
For values of the transverse momentum Qt greater than roughly half the pair 
mass Q, Qt > (5/2, we show that the differential cross section is dominated by 
subprocesses initiated by incident gluons, provided that the polarized gluon 
density is not too small. Massive lepton-pair differential cross sections should 
be a good source of independent constraints on the polarized gluon density, 
free from the experimental and theoretical complications of photon isolation 
that beset studies of prompt photon production. We provide predictions for 
the spin-averaged and spin-dependent differential cross sections as a function 
of Qt at energies relevant for the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven, and we compare these with predictions for real prompt photon 
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production. 
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I. INTRODUCTION AND MOTIVATION 



Both massive lepton-pair production, /ii + /i2 ^ 7* + X; 7* — > II, and prompt real 
photon production, hi + h2 ^ ■J + X are valuable probes of short-distance behavior in 
hadron reactions. The two reactions supply critical information on parton momentum den- 
sities, in addition to the opportunities they offer for tests of perturbative quantum chro- 
modynamics (QCD). Spin-averaged parton momentum densities may be extracted from 
spin-averaged nucleon-nucleon reactions, and spin-dependent parton momentum densities 
from spin-dependent nucleon-nucleon reactions. An ambitious experimental program of 
measurements of spin-dependence in polarized proton-proton reactions will begin soon at 
Brookhaven's Relativistic Heavy Ion Collider (RHIC) with kinematic coverage extending 
well into the regions of phase space in which perturbative quantum chromodynamics should 
yield reliable predictions. 

Massive lepton-pair production, commonly referred to as the Drell-Yan process pro- 
vided early confirmation of three colors and of the size of next-to-leading contributions to 
the cross section differential in the pair mass Q. The mass and longitudinal momentum (or 
rapidity) dependences of the cross section (integrated over the transverse momentum Qt 
of the pair) serve as laboratory for measurement of the antiquark momentum density, com- 
plementary to deep-inelastic lepton scattering from which one gains information of the sum 
of the quark and antiquark densities. Inclusive prompt real photon production is a source 
of essential information on the gluon momentum density. At lowest order in perturbation 
theory, the reaction is dominated at large values of the transverse momentum pt of the 
produced photon by the "Compton" subprocess, q + g + q- This dominance is preserved 
at higher orders, indicating that the experimental inclusive cross section differential in pt 
may be used to determine the density of gluons in the initial hadrons [@-|^. 

In two previous papers 0, we addressed the production of massive lepton-pairs as a func- 
tion of the transverse momentum Qt of the pair in unpolarized nucleon-nucleon reactions, 
/ii + /i2 7* + X, in the region where Qt is greater than roughly half of the mass of the 
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pair, Qt > Q/2. We demonstrated that the differential cross section in this region is dom- 
inated by subprocesses initiated by incident gluons. Correspondingly, massive lepton-pair 
differential cross sections in unpolarized nucleon-nucleon reactions are a valuable, heretofore 
overlooked, independent source of constraints on the spin-averaged gluon density. 

Turning to longitudinally polarized proton-proton collisions in this paper, we study the 
potential advantages that the Drell-Yan process may offer for the determination of the 
spin-dependence of the gluon density. To be sure, the cross section for massive lepton-pair 
production is smaller than it is for prompt photon production. However, just as in the 
unpolarized case, massive lepton pair production is cleaner theoretically since long-range 
fragmentation contributions are absent as are the experimental and theoretical complications 
associated with isolation of the real photon. Moreover, the dynamics of spin-dependence in 
hard-scattering processes is a sufficiently complex topic, and its understanding at an early 
stage in its development, that several defensible approaches for extracting polarized parton 
densities deserve to be pursued with the expectation that consistent results must emerge. 

There are notable similarities and differences in the theoretical analyses of massive lepton- 
pair production and prompt real photon production. At first-order in the strong coupling 
strength, a^, the Compton subprocess and the annihilation subprocess q + q ^ 'J + g supply 
the transverse momentum of the directly produced prompt photons. Identical subprocesses, 
with the real 7 replaced by a virtual 7*, are responsible at 0{as) for the transverse mo- 
mentum of massive lepton-pairs. An important distinction, however, is that fragmentation 
subprocesses play a very important role in prompt real photon production at collider en- 
ergies. In these long-distance fragmentation subprocesses, the photon emerges from the 
fragmentation of a final parton, e.g., q + g q + g, followed by g 7 -|- X. The neces- 
sity to invoke phenomenological fragmentation functions and the infrared ambiguity ^ of 
the isolated cross section in next-to-leading order raise questions about the extent to which 
isolated prompt photon data may be used for fully quantitative determinations of the gluon 
density. It is desirable to investigate other physical processes for extraction of the gluon 
density that are free from these systematic uncertainties. Fortunately, no isolation would 
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seem necessary in the case of virtual photon production (and subsequent decay into a pair of 
muons) in typical collider or fixed target experiments. Muons are observed only after they 
have penetrated a substantial hadron absorber. Thus, any hadrons within a typical cone 
about the direction of the 7* will have been stopped, and the massive lepton-pair signal will 
be entirely inclusive. 

Another significant distinction between massive lepton-pair production and prompt real 
photon production is that interest in /ii + /i2 — > 7* + X has been drawn most often to the 
domain in which the pair mass Q is relatively large, justifying a perturbative treatment based 
on a small value of as{Q) and the neglect of inverse-power high-twist contributions (except 
near the edges of phase space). The focus in prompt real photon production is directed to 
the region of large values of pt where as(pr) is small. Interest in the transverse momentum 
Qt dependence of the massive lepton-pair production cross section has tended to be limited 
to small values of Qt where the cross section is largest. Fixed-order perturbation theory |^ 
is applicable for large Qt, but it is inadequate at small Qt, and all-orders resummation 
methods |§-|T3| have been developed to address the region Qt « Q- 

As long as Qt is large, the perturbative requirement of small q;s(Qt) can be satisfied 
without a large value of Q. We therefore explore and advocate the potential advantages of 
studies of (Pa/dQdQT as a function of Qt for modest values of Q, <5 ~ 2 to 3 GeV, below 
the range of the traditional Drell-Yan region. There are various backgrounds with which to 
contend at small Q such as the contributions to the event rate from prompt decays of heavy 
flavors, e.g., hi + h2 — ^ c + c + X; c ^ / + X. These heavy flavor contributions may be 



estimated by direct computation and/or bounded through experimental measurement 
of the like-sign-lepton distributions. 

In Sec. II, we review perturbative QCD calculations of the transverse momentum dis- 
tribution for massive lepton-pair production in the case in which the initial nucleon spins 
are polarized as well as in the spin-average case. In Sec. Ill, we present next-to-leading 
order predictions for the transverse momentum dependence of the cross sections for mas- 
sive lepton-pair and real prompt photon production in unpolarized proton-proton collisions 



at energies typical of the RHIC collider. Predictions for spin dependence are provided in 
Sec. IV. Our conclusions are summarized in Sec. V. 



II. MASSIVE LEPTON PAIR PRODUCTION AND PROMPT PHOTON 
PRODUCTION AT NEXT-TO-LEADING ORDER 

In inclusive hadron interactions at collider energies, + — ^ 7*+-^^ with 7* II, lepton 
pair production proceeds through partonic hard-scattering processes involving initial-state 
light quarks q and gluons g. In lowest-order QCD, at 0{a^), the only partonic subprocess is 
g + g — > 7*. At 0{as), both g + g — > 7* + (? and g + (? — > 7* + g participate, with the recoil of 
the final parton balancing the transverse momentum of the lepton-pair. These processes are 
shown in Figs. 1(a) and 2(a). Calculations of the cross section at order 0{al) involve virtual 
loop corrections to these 0{as) subprocesses (Figs. 1(b) and 2(b)) as well as contributions 
from a wide range of 2 — > 3 parton subprocesses (of which some examples are shown in 
Figs. 1(c) and 2(c)). 

The physical cross section is obtained through the factorization theorem, 

^ = E / dx,dx,fi{x„ ^))fi{x,, /xj)^(., g, Qt, y; /xj). (1) 

It depends on the hadronic center-of-mass energy S and on the mass Q, the transverse 
momentum Qt, and the rapidity y of the virtual photon; /x/ is the factorization scale of the 
scattering process. The usual Mandelstam invariants in the partonic system are defined by 
•5 — (Pi +^2)^, t — {pi — P-y*)"^-! and u — {p2 — P-y*)"^, where pi and p2 are the momenta 
of the initial state partons and p^* is the momentum of the virtual photon. The indices 
ij e {qq, qg} denote the initial parton channels whose contributions are added incoherently 
to yield the total physical cross section. Functions fl{x,ii) denote the usual spin-averaged 
parton distribution functions. 

The partonic cross section a]j (s, Q, Qt, y, jJ^f) is obtained commonly from fixed-order 
QCD calculations through 
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(fa'^* (fa'^*'^"^ c^2aT*'(^) (prr^^^"^ 

dtdu dtdu dtdu dtdu 

The tree, virtual loop, and real emission contributions are labeled (a), (b), and (c) as are 
the corresponding diagrams in Figs. 1 and 2. The parameter /i is the renormalization scale. 



It is set equal to the factorization scale /x/ = yQ^ + Qj' throughout this paper. 

The cross section for hi + h2 —>■ II + X , differential in the invariant mass of the lepton 
pair as well as its transverse momentum and rapidity, is obtained from Eq. ([l|) by the 
relation 

-I ^^-\^^^^'^^^S,Q,Qr,y), (3) 



dQ'^dQj.dy \3ttQ'^ I dQj,dy 

where = [pi +PiY, and pi,Pi are the four- momenta of the two final leptons. The Drell-Yan 
factor aem/ (SttQ^) is included in all numerical results presented in this paper. 

While the full next-to-leading order QCD calculation exists for massive lepton-pair pro- 
duction in the case of unpolarized initial nucleons, only a partial calculation is available in the 
polarized case [|1^]. Correspondingly, we present spin-averaged differential cross sections at 
next-to-leading order, but we calculate spin asymmetries at leading order. Spin asymmetries 
are obtained by dividing the spin-dependent differential cross section by its spin-averaged 
counterpart. For prompt photon production, comparisons of asymmetries computed at next- 
to- leading order with those at leading order show only modest differences whereas the 
cross sections themselves are affected more significantly. Given the similarity of prompt 
photon production and massive lepton-pair production in the region of Qt of interest to 
us 0, we expect that the leading-order asymmetries will serve as a useful guide for massive 
lepton-pair production. 

Rewriting Eq. (Q) and integrating over an interval in Q'^, we calculate the spin-averaged 
differential cross section Ed^a1^_^y^^/ dp^ as 



dp^ ~ Sn^Sj-JQ^ . J _ hAxi,f,j)h^{x2,fij)s . (4) 



In Eq. (Q), Q'^a^ and (5^j„ are the chosen upper and lower limits of integration for Q"^, and 



^mm _ j^^^ _ ^y(^i _ jf.^y '-£\^Q value of X2 is determined from X2 = {xiX2 — t)/{xi — xi), 
with 



^1 = ^-g^ = ^eV4 + 4r, (5) 



and 



X2 = = \e-y^^^. (6) 

We use Pi and P2 to denote the four-vector momenta of the incident nucleons; S = (Pi+P2)^- 
The invariants in the hadronic system, T = (Pi — p-y^ and U = (P2 — Pj*)"^, are related to 
the partonic invariants by 

(t-Q^)=xi{T- = -X1X2S, (7) 

and 

(n - Q2) = a;2(f/ - Q') = -xaxi^. (8) 
The scaled variables xt and r are 



When the initial nucleons are polarized longitudinally, we can compute the difference of 
cross sections 

Act = a++ -(T+_, (9) 

where +, — denote the helicities of the incident nucleons. 
In analogy to Eq. (|), we find 

The functions Afl{x,^) denote the spin-dependent parton distribution functions, defined 
by 
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(11) 



fl^{x,fif) is the distribution of partons of type i with positive (+) or negative (— ) hehcity 
in hadron h. Likewise, the polarized partonic cross section Aa"'' is defined by 



Aa^*=a^*(+,+)-a^*(+,-), 



(12) 



with +, — denoting the hehcities of the incoming partons. 

The hard subprocess cross sections in leading order for the unpolarized and polarized 
cases are 



'qq 



dt 



' dt Nc s 



u t_ 2g^(g^ - M-t) 

t u ut 



(13) 



and 



da. 
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dt 



2 TTttem (Xg 



Q — = P 

dt ^ Nc s 



s t 2Q'^u 
t s st 



2u + s 2u + t 



t 



(14) 



(15) 



Our results on the longitudinal spin dependence are expressed in terms of the two-spin 
longitudinal asymmetry All, defined by 



A 



a • 



LL 



where +, — denote the helicities of the incoming protons. 



(16) 



III. UNPOLARIZED CROSS SECTIONS 

We turn in this Section to explicit evaluations of the differential cross sections as func- 
tions of Qt at collider energies. We work in the MS renormalization scheme and set the 
renormalization and factorization scales equal. We employ the MRST set of spin-averaged 
parton densities [IT^ and a two-loop expression for the strong coupling strength Q;s(/i), with 



five flavors and appropriate threshold behavior at /i = mj,; A*^^) = 300 MeV. The strong 



coupling strength ag is evaluated at a hard scale fj, — yjQ'^ + Qt- We present results for 
three values of the ccnter-of-mass energy, = 50, 200, and 500 GeV. 

For — 200 GeV, we present the invariant inclusive cross section Ed^a/dp^ as a 
function of Qt in Fig. 3. Shown in this figure are the qq and qg perturbative contributions 
to the cross section at leading order and at next-to-leading order. We average the invariant 
inclusive cross section over the rapidity range -1.0 < y < 1.0 and over the mass interval 5 
< Q < 6 GeV. For Qt < 1.5 GeV, the qq contribution exceeds that of qg channel. However, 
for values of Qt > 1.5 GeV, the qg contribution becomes increasingly important. As shown 
in Fig. 4(a), the qg contribution accounts for about 80 % of the rate once Qt — Q- The 
results in Fig. 4(a) also demonstrate that subprocesses other than those initiated by the qq 
and qg initial channels contribute negligibly. 

In Fig. 4(b), we display the fractional contributions to the cross section as a function 
of Qt for a larger value of Q: 11 < Q < 12 GeV. In this case, the fraction of the rate 
attributable to qg initiated subprocesses again increases with Qt- It becomes 80 % for 
Qtc^Q. 

For the calculations reported in Figs. 3 and 4(a,b), we chose values of Q in the traditional 
range for studies of massive lepton-pair production, viz., above the interval of the J/ip and 
ip' states and either below or above the interval of the T's. 

For Fig. 4(c), we select the interval 2.0 < Q < 3.0 GeV. In this region, one would be 
inclined to doubt the reliability of leading- twist perturbative descriptions of the cross section 
da/dQ, integrated oyer all Qt- However for values of Qt that are large enough, a perturbative 
description of the Qt dependence of d'^a/dQdQT ought to be justified. The results presented 
in Fig. 4(c) demonstrate that, as at higher masses, the qg incident subprocesses dominate 
the cross section for Qt — Q- 

The calculations presented in Figs. 4 show convincingly that data on the transverse 
momentum dependence of the cross section for massive lepton-pair production at RHIC 
collider energies should be a valuable independent source of information on the spin-averaged 
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gluon density. 

In Fig. 5, we provide next-to- leading order predictions of the differential cross section as 
a function of Qt for three values of the center-of-mass energy and two intervals of mass Q. 
Taking Ed^a/dp^ = 10-=^pb/GeV^ as the minimum accessible cross section, we may use the 
curves in Fig. 5 to establish that the massive lepton-pair cross section may be measured to 
Qt = 7.5, 14, and 18.5 GeV at \/5 = 50, 200, and 500 GeV, respectively, when 2 <Q < 
3 GeV, and to Qt — 6, 11.5, and 15 GeV when 5 < Q < 6 GeV. In terms of reach in the 
fractional momentum Xgiuon carried by the gluon, these values of Qt may be converted to 
Xgiuon ^ XT = 2Qt/^ = 0.3, 0.14, and 0.075 at ^/S = 50, 200, and 500 GeV when 2 
< g < 3 GeV, and to Xgiuon - 0.24, 0.115, and 0.06 when 5 <Q <6 GeV. On the face of 
it, the smallest value of \/S provides the greatest reach in Xgiuon- However, the reliability 
of fixed-order perturbative QCD as well as dominance of the qg subprocess improve with 
greater Qt- The maximum value Qt — 7.5 GeV attainable at = 50 GeV argues for a 
larger y/S. 

It is instructive to compare our results with those expected for prompt real photon pro- 
duction. In Fig. 6, we present the predicted differential cross section for prompt photon 
production for three center-of-mass energies. We display the result with full fragmenta- 
tion taken into consideration (upper line) and with no fragmentation contributions included 
(lower line). Comparing the magnitudes of the prompt photon and massive lepton pair 
production cross sections in Figs. 5 and 6, we note that the inclusive prompt photon cross 
section is a factor of 1000 to 4000 greater than the massive lepton-pair cross section inte- 
grated over the mass interval 2.0 < Q < 3.0 GeV, depending on the value of Qt- This 
factor is attributable in large measure to the factor ctem/ (37r(5^) associated with the decay 
of the virtual photon to /x^/i.^. Again taking Ed^a/dp^ = 10^'^pb/GeV^ as the minimum 
accessible cross section, we may use the curves in Fig. 6 to establish that the real photon 
cross section may be measured to pt — 14, 33, and 52 GeV at = 50, 200, and 500 GeV, 
respectively. The corresponding reach in xt = 2pT/y/S = 0.56, 0.33, and 0.21 at = 50, 
200, and 500 GeV is two to three times that of the massive lepton-pair case. 
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The breakdown of the real photon direct cross section at \^ = 200 GeV into its qq and 
qg components is presented in Fig. 7. As may be appreciated from a comparison of Figs. 4 
and 7, dominance of the qg contribution in the massive lepton-pair case is as strong as in the 
prompt photon case. The significantly smaller cross section in the case of massive lepton- 
pair production means that the reach in Xgiuon is restricted to about a factor of two to three 
less, depending on and Q, than that potentially accessible with prompt photons in the 
same sample of data. Nevertheless, it is valuable to be able to investigate the gluon density 
with a process that has reduced experimental and theoretical systematic uncertainties from 
those of the prompt photon case. 

In our previous papers we compared our spin-averaged cross sections with available 
fixed-target and collider data on massive lepton-pair production at large values of Qt, and 
we were able to establish that fixed-order perturbative calculations, without resummation, 
should be reliable for > Q/2. The region of small Qt and the matching region of 
intermediate Qt are complicated by some level of phenomenological ambiguity. Within 
the resummation approach, phenomenological non-perturbative functions play a key role 
in fixing the shape of the Qt spectrum at very small Qt, and matching methods in the 
intermediate region are hardly unique. For the goals we have in mind, it would appear best 
to restrict attention to the region Qt > Q/2. 



IV. PREDICTIONS FOR SPIN DEPENDENCE 

Given theoretical expressions derived in Sec. II that relate the spin-dependent cross 
section at the hadron level to spin-dependent partonic hard-scattering matrix elements and 
polarized parton densities, we must adopt models for spin-dependent parton densities in 
order to obtain illustrative numerical expectations. For the spin-dependent densities that 
we need, we use the three different parametrizations suggested by Gehrmann and Stirling 



(GS) ||T8[. We have verified that the positivity requirement Afl{x, fif)/ fl{x, fif] 
satisfied. 



< 1 is 
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The current deep inelastic scattering data do not constrain the polarized gluon density 
tightly, and most groups present more than one plausible parametrization. Gehrmann and 
Stirling present three such parametrizations, labelled GSA, GSB, and GSC. In the GSA 



and GSB sets, AG(a;, /lo) is positive for all x, whereas in the GSC set AG{x, fio) changes 
sign. After evolution to fi'j = 100 GeV^, AG{x,fif) remains positive for essentially all x in 
all three sets, but its magnitude is small in the GSB and GSC sets. 

In this Section, we present two-spin longitudinal asymmetries for massive lepton-pair 
production as a function of transverse momentum. Results are displayed for pp collisions at 
the center-of-mass energies = 50, 200, and 500 GeV typical of the Brookhaven RHIC 
collider. 

In Figs. 8(a-c), we present the two-spin longitudinal asymmetries. All, as a function of 
Qt- As noted earlier, these asymmetries are computed in leading-order. More specifically, 
we use leading-order partonic subprocess cross sections a and Aa with next-to-leading order 
spin-averaged and spin-dependent parton densities and a two- loop expression for a^. The 
choice of a leading-order expression for Aa is required because the full next-to-leading order 
derivation of Aa has not been completed for massive lepton-pair production. Experience 
with prompt photon production indicates that the leading-order and next-to-leading order 
results for the asymmetry are similar so long as both are dominated by the qg subprocess. 
Results are shown for three choices of the polarized gluon density. The asymmetry becomes 
sizable for large enough Qt for the GSA and GSB parton sets but not in the GSC case. 
Comparing the three figures, we note that All is nearly independent of the pair mass Q as 
long as Qt is not too small. This feature should be helpful for the accumulation of statistics; 
small bin- widths in mass are not necessary, but the J/ ip and T resonance regions should be 
excluded. 

As noted above the qg subprocess dominates the spin- averaged cross section. It is in- 
teresting and important to inquire whether this dominance persists in the spin-dependent 
situation. In Figs. 9 and 10, we compare the contribution to the asymmetry from the polar- 
ized qg subprocess with the complete answer for all three sets of parton densities. The qg 
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contribution is more positive than the full answer for values of Qt that are not too small; 
the full answer is reduced by the negative contribution from the qq subprocess for which 
the parton-level asymmetry cill = — 1- At small Qt-, the net asymmetry may be driven 



negative by the qq contribution, and based on our experience with other calculations |]T9 
from processes such as gg that contribute in next-to-leading order. For the GSA and GSB 
sets, we see that once it becomes sizable (e.g., 5% or more), the total asymmetry from all 
subprocesses is dominated by the large contribution from the qg subprocess. 

As a general rule in studies of polarization phenomena, many subprocesses can contribute 
small and conflicting asymmetries. Asymmetries are readily interpretable only in situations 
where the basic dynamics is dominated by one major subprocess and the overall asymmetry is 
sufficiently large. In the case of massive lepton-pair production that is the topic of this paper, 
when the overall asymmetry An itself is small, the contribution from the qg subprocess 
cannot be said to dominate the answer. However, if a large asymmetry is measured, similar 
to that expected in the GSA case at the larger values of Qt-, Figs. 9 and 10 show that the 
answer is dominated by the qg contribution, and data will serve to constrain AG{x, fif). If 
AG(x,/i/) is small and a small asymmetry is measured, such as for the GSC parton set, or 
at small Qt for all parton sets, one will not be able to conclude which of the subprocesses is 
principally responsible, and no information could be adduced about AG(x, /i/), except that 
it is small. 

In Figs. 10 (a) and (b), we examine the energy dependence of our predictions for two dif- 
ferent intervals of mass Q. For Qt not too small, we observe that Alt in massive lepton pair 
production is well described by a scaling function of xt = 2Qt/ VS, All{^/S, Qt) — h^* (xt)- 
In our discussion of the spin-averaged cross sections, we took Ed^cr/dp^ = 10^'^pb/GeV^ as 
the minimum accessible cross. Combining the results in Fig. 5 with those in Fig. 10, we see 
that longitudinal asymmetries An = 20%, 7.5%, and 3% are predicted at this level of cross 
section at = 50, 200, and 500 GeV when 2 < Q < 3 GeV, and All = 11%, 5%, and 2% 
when 5 < Q < 6 GeV. For a given value of Qt, smaller values of result in greater 
asymmetries because AG{x)/G{x) grows with x. 
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The predicted cross sections in Fig. 5 and the predicted asymmetries in Fig. 10 should 
make it possible to optimize the choice of center-of-mass energy at which measurements 
might be carried out. At \fS = 500 GeV, asymmetries are not appreciable in the interval 
of Qt in which event rates are appreciable. At the other extreme, the choice of = 50 
GeV does not allow a sufficient range in Qt- Accelerator physics considerations favor higher 
energies since the instantaneous luminosity increases with v^- Investigations in the energy 
interval = 150 to 200 GeV would seem preferred. 

In Fig. 11, we display predictions for All in prompt real photon production for three 
values of the center-of-mass energy. These calculations are done at next-to-leading order in 
QCD. Dominance of the qg contribution is again evident as long as All is not too small. 
So long as Qt > Q, we note that the asymmetry in massive lepton-pair production is 
about the same size as that in prompt real photon production, as might be expected from 
the strong similarity of the production dynamics in the two cases. As in massive lepton- 
pair production. All in prompt photon production is well described by a scaling function 
of XT = 2pt/VS, ALLi^: Pt) ~ hy{xT)- For Ed?a/dp^ = 10 ^pb/GeV^, we predict 
longitudinal asymmetries All — 31%, 17%, and 10% in real prompt photon production at 
VS = 50, 200, and 500 GeV. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we focus on the Qt distribution for — > 7*-|-X. We present and discuss 
calculations carried out in QCD at RHIC coUider energies. We show that the differential 
cross section in the region Qt > Q/2 is dominated by subprocesses initiated by incident 
gluons. Dominance of the qg contribution in the massive lepton-pair case is as strong as 
in the prompt photon case, p + p ^ ^ + X . As our calculations demonstrate, the Qt 
distribution of massive lepton pair production offers a valuable additional method for direct 
measurement of the gluon momentum distribution. The method is similar in principle to the 
approach based on prompt photon production, but it avoids the experimental and theoretical 
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complications of photon isolation that beset studies of prompt photon production. 

As long Qt is large, the perturbative requirement of small ag^Qr) can be satisfied without 
a large value of Q. We therefore explore and advocate the potential advantages of studies of 
(fa /dQdQx as a function of Qt for modest values of Q, Q ~ 2 GeV, below the traditional 
Drell-Yan region. 

For the goals we have in mind, it would appear best to restrict attention to the region in 
Qt above the value at which the resummed result falls below the fixed-order perturbative 
expectation. A rough rule-of-thumb based on our calculations is Qt > Q/2. Uncertainties 
associated with resummation make it impossible to use data on the Qt distribution at small 
Qt to extract precise information on parton densities. 

In this paper we also present a calculation of the longitudinal spin-dependence of mas- 
sive lepton-pair production at large values of transverse momentum. We provide polarization 
asymmetries as functions of transverse momenta that may be useful for estimating the fea- 
sibility of measurements of spin-dependent cross sections in future experiments at RHIC 
collider energies. The Compton subprocess dominates the dynamics in longitudinally po- 
larized proton- proton reactions as long as the polarized gluon density AG{x,iif) is not too 
small. As a result, two-spin measurements of inclusive prompt photon production in polar- 
ized pp scattering should constrain the size, sign, and Bjorken x dependence of AG{x,iJ,f). 
Significant values of An (i.e., greater than 5 %) may be expected for xt — 2Qt/VS > 0.10 
if the polarized gluon density AG{x,iJ,f) is as large as that in the GSA set of polarized parton 
densities. If so, the data could be used to determine the polarization of the gluon density 
in the nucleon. On the other hand, for small AG{x,iif), dominance of the qg subprocess is 
lost, and AG'(x, is inaccessible. 
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(a) 




(b) 





FIG. 1. (a)Lowest-order Feynman diagrams for the direct process q+q — 7+5. (b) Examples of 
virtual gluon loop diagrams, (c) Examples of next-to-leading order three-body final-state diagrams. 
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FIG. 2. As in Fig. 1, but for the subprocesses initiated by the q + g initial state. 
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pp yX at VS = 200 GeV 




5 10 15 20 25 30 

Qt (GeV) 



FIG. 3. Lowest order (dashed lines) and next-to-leading order (solid lines) perturbative calcu- 
lations of the invariant inclusive cross section Ed^a/dp^ as a function of Qt for pp 7* 
200 GeV, in the MS scheme. Contributions from the qg and qq channels are shown separately. The 
results are averaged over the rapidity interval - 1.0 < y < 1.0 and over the interval 5.0 < Q < 6.0 
GeV. 
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pp yX at VS = 200 GeV 
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FIG. 4. Contributions from the various partonic subprocesses to the invariant inclusive cross 
section Ed^a/dp^ as a function of Qt for pp — ^ j*X at = 200 GeV. The cross section is 
averaged over the rapidity interval -1.0 < y < 1.0 and over the intervals (a) 5.0 < Q < 6.0 GeV, 
(b) 11.0 < Q < 12.0 GeV, and (c) 2.0 < Q < 3.0 GeV. The contributions are labeled by qg (solid), 
qq (dashed), gg (dotted), qq2 non-factorizable parts (dot-dashed), and qq (wide dots). 
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FIG. 5. Differential cross sections Ed^a/dp^ as a function of Qt for for pp ^ ^* + X &t s/S = 
50, 200, and 500 GeV, averaged over the rapidity interval -1.0 < y < 1.0 and the mass intervals 
(a) 2.0 < g < 3.0 GeV, and (b) 5.0 < Q < 6.0 GeV. 
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FIG. 6. Differential cross sections Ed^a/dp^ at rapidity y = as a function of pr for real 
photon production pp at VS = 50, 200, and 500 GeV for two cases: no fragmentation terms 

included (lower) and the inclusive case with full fragmentation included (upper). 
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FIG. 7. Fractional contributions from the next-to-leading order qq and the qg channels to the 
differential cross sections Ed^a/dp^ as a function of px for real photon production pp — 7X at 
a/5 = 200 GeV and rapidity y = 0. No fragmentation terms are included. 
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pp yX at VS = 200 GeV 
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FIG. 8. Computed longitudinal asymmetry All as a function of Qt for pp —>■ j*X at s/S = 200 
GeV. The asymmetry is averaged over the rapidity interval -1.0 < y < 1.0 and over the intervals 
(a) 2.0 < Q < 3.0 GeV, (b) 5.0 < Q < 6.0 GeV, and (c) 11.0 <Q < 12.0 GeV, for three choices 
of the polarized parton densities, GSA, GSB, and GSC. 
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FIG. 9. Comparison of the contribution of the qg subprocess (upper curve) to the longitudinal 
asymmetry with the total (lower curve) as a function of Qt for -p-p 7* X at ^/S = 200 
GeV. The asymmetry is averaged over the rapidity interval -1.0 < y < 1.0 and over the interval 
5.0 < Q < 6.0 GeV. Results are shown for three sets of spin-dependent parton densities. 
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FIG. 10. Computed longitudinal asymmetry as a function of Qt for for pp ^ ^* + X at 
VS = 50, 200, and 500 GeV, averaged over the rapidity interval -1.0 < y < 1.0 and the mass 
intervals (a) 5.0 < Q < 6.0 GeV and (b) 2.0 < Q < 3.0 GeV. Shown are both the complete answer 
at leading-order and the contribution from the qg subprocess. The GSA set of polarized parton 
densities is used. 
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FIG. 11. Computed longitudinal asymmetry A^,/, as a function of for real photon production 
pp — ^ 7 + X at \fS = 50, 200, and 500 GeV and rapidity y = 0. Shown are both the complete 
answer at leading-order and the contribution from the qg subprocess. The GSA set of polarized 
parton densities is used. 
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